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ABSTRACT: Reliability-based design optimization (RBDO) has gained much attention in the past fifteen 



Oscars as a way of introducing robustness in the process of designing structures and systems in an optimal manner. 
Indeed classical optimization (e.g. minimize some cost under mechanical constraints) usually leads to solutions 
that lie at the boundary of the admissible domain, and that are consequently rather sensitive to uncertainty in the 
^design parameters. In contrast, RBDO aims at designing the system in a robust way by minimizing some cost 
function under reliability constraints. Thus RBDO methods have to mix optimization algorithms together with 
rTjl'eliability calculations. The classical approach known as "double-loop" consists in nesting the computation of 
Sthe failure probability with respect to the current design within the optimization loop. It is not applicable to 
1^ industrial models (e.g. finite element models) due to the associated computational burden. In contrast, methods 
+^based on the approximation of the reliability (e.g. FORM) may not be applicable to real-world problems. In 
-^this context, an original method has been developed that tries to circumvent the abovementioned drawbacks of 
I— the existing approaches. It is based on the adaptive construction of a meta-model for the expensive-to-evaluate 
^mechanical model, and on the subset simulation technique for the efficient and accurate computation of the 
j>failure probability and its sensitivities with respect to the design variables. The proposed methodology is briefly 
0>described in this paper before it is applied to the reliability-based design of an imperfect submarine pressure 
fThuU. 
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INTRODUCTION 



;:*In order to ensure that their design is safe enough, 
• (mechanical engineers and other stake-holders usu- 
^ally cope with uncertainty using partial safety factors. 
^Such factors are usually prescribed in specific stan- 
dards depending on the field of application and are 
calibrated with respect to both experience and scat- 
ter of environmental, material and structural proper- 
ties. However, they might either introduce an exces- 
sive degree of conservatism or on the contrary, they 
might lack exhaustivity about the servicing condi- 
tions of the system to be designed. This is the rea- 
son why Reliability-Based Design Optimization has 
gained much attention in the past decades as it allows 
one to consider explicitly uncertainty in the design op- 
timization procedure. 

RBDO aims at designing the system in a robust 
manner by minimizing some cost function c with re- 
spect to a set of characteristic design parameters 
(e.g. means) under reliability constraints instead of 



the original deterministic constraints. The reliability- 
based optimal design 6* of interest is thus defined as 
the following minimizer: 



9* = axs min c (9) 



fi{9)<0,i = l,...,nc 
F{gi(X(9))<0)<P^l,l 



(1) 



where Vg is the admissible design space, {fi,i = 
1, ...,nc} are classical deterministic constraints 
that either do not put the system's safety at 
stake or involves partial safety factors in order 
to cope with some unconsidered failure scenarii, 
{P (gi {X (0)) < 0) < P^;, ; = 1, rip} are the so-called 
reliability constraints involving the vector of uncer- 
tain variables X and the maximum failure proba- 
bilities that are tolerated by the code of practice 



,np}. Note that, in this paper, X is dis- 
tributed according to its joint probability density func- 
tion fx (•, 9) which is parameterized by the design pa- 
rameters 9. 



In essence, RBDO methods have to mix optimiza- 
tion algorithms together with reliability calculations. 
The classical approach known as ''double loop" con- 
sists in nesting the computation of the probability 
of failure with respect to the current design within 
the optimization loop. Crude or advanced simulation 
methods are not applicable to industrial models {e.g. 
finite element models) due to the associated compu- 
tational burden. In contrast, methods based on the ap- 
proximation of the reliability {e.g. FORM) may not 
be applicable to real- world problems. The interested 
reader may refer to the book by jTsompanakis et aL 
([2008) which provides a complete review of existing 
approaches. 

In order to try to circumvent the aforementioned 
drawbacks of existing approaches, an original method 
has been developed. It is based on both the adaptive 
construction of a meta-model for the original time- 
consuming mechanical model, and the use of the sub- 
set simulation technique for the efficient and accurate 
computation of the failure probability and its associ- 
ated sensitivities w.r.t. the design variables. The strat- 
egy is briefly described in this paper before its appli- 
cation to the design of imperfect submarine pressure 
hulls. The interested reader may find some more de- 
tailed explanations in a forthcoming paper (Dubourg 
et al.[ 2011). The case addressed in this paper aims 
at finding the optimal dimensions of a single bay of 
a ring-stiffened cylinder representative of a subma- 
rine pressure hull under external hydrostatic pressure. 
It accounts for uncertainties in the material proper- 
ties and in the amplitudes of its imperfections. These 
imperfections are supposed distributed over the two 
most critical buckling modes in a first simplified ap- 
proach (overall and interframe buckling modes). In- 
deed, the work presented in this paper partly relies 
on closed-form and semi-analytical formulas that are 
commonly used for preliminary design in the subma- 
rine industry {e.g. the ,BS5500) . The results provided 
in this paper should provide a first sound basis for 
a future application involving a higher fidelity non- 
linear shell finite element model. 

2 THE KRIGING SURROGATE 

In computer experiments, surrogate modelling ad- 
dresses the problem of replacing a time-consuming 
mathematical model M called a simulator by an em- 
ulator that is much faster to evaluate. The emulator 
has the same input space Vx c and output space 
c M as the simulator. For the sake of simplicity 
only scalar-output simulators are considered here. It 
is important to note that both x and y are determinis- 
tic in this context. The emulator is usually built from 
a set of observations named a design of experiments 
(DOE) {yi =M{xi), i = \y .m }- 

Kriging (Santner et al. 2003 [ ) is one particular em- 
ulator that is able to give a probabilistic response 
Y {x) whose variance (scatter) depends on the quan- 



tity of available knowledge. In other words, the un- 
certainty in the prediction is purely epistemic and thus 
reducible as shown in the next sections. 

In essence, kriging assumes that the output y = 
M {x) is a sample path of an underlying Gaussian 
stochastic process Y [x) that would read as follows: 



Y{x)^ f{x)^ f5 + Z{x) 

where: 



(2) 



• / [x]^ (3 denotes the mean of the GP which cor- 
responds to a classical linear regression model 
on a given functional basis {/j, i = ... ,p} e 

• Z (x) denotes the stochastic part of the GP which 
is modelled as a zero mean, constant variance o-y , 
stationary Gaussian process with a given sym- 
metric positive definite autocorrelation model. It 
is fully defined by its autocovariance function 
which reads {{x, x') eV^ y< Vx): 

Cyy{x,x') =a'^R{\x-x'\,£) (3) 

where £ is a vector of parameters defining R. 

The most widely used class of autocorrelation func- 
tions is the anisotropic squared exponential model: 




(4) 



The best linear unbiased estimation of the GP Y at 



point X is shown ( Santner et al. 2003 ) to be the Gaus- 
sian random variate Y {x) ~ N {y {x) , {x)^ with 
moments: 



Y{x) = f{xf f3 + r{xfR-^ (y -F^ 



(5) 



where we have introduced r, R and F such that: 
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(6) 



n = R{\x-Xi\,£), i^l, ... ,m 
Rij = R(\xi-Xj\,£) , i = l, ... ,m, j ^1, 



Fi. 



fi {•^i 



, m 



(7) 

,m (8) 

(9) 



At this stage one can easily prove that Y {xi) = M (xi) 
and crj; (a;j) = for i = 1, ... , m, thus meaning the krig- 
ing surrogate is an exact interpolator. 

Given a choice for the regression and correlation 
models, the optimal set of parameters /3*, £* and 
ai* can then be inferred using the maximum likeli- 
hood principle applied to the unique sparse obser- 
vation of the GP sample path grouped in the vector 



y = {M {xi 



This optimization problem 



can be solved analytically for both f3* and a'l* assum- 
ing £* is known. Thus, the problem is solved in two 
steps: the maximum likelihood estimation of £* is first 
solved by a global optimization algorithm which in 
turns allows one to evaluate the optimal (3* and a'i*. 



3 THE PROPOSED ADAPTIVE REFINEMENT 
TECHNIQUE 

Various refinement techniques have been proposed in 
the kriging-related literature, especially in structural 
reliability ( |Bichon et alj [20U8| |Lee and Jungj [2008] ). 
They are all based on the genuine probabilistic na- 
ture of the kriging prediction. The refinement tech- 
nique consists in finding the set of points that should 
be added to the DOE in order to improve the predic- 
tion and reduce its associated epistemic uncertainty in 
a specific region of interest. 

In reliability and reliability-based design optimiza- 
tion, the region of interest in the space of random vari- 
ables X is the limit-state surface g = - we drop the 
subscript / in this section for the sake of clarity. Pro- 
vided an initial kriging prediction of the performance 
function g, the probable vicinity of the limit-state sur- 
face g = can be defined in terms of the following 
margin of uncertainty: 



|a; : —ka^ [x) < G (x) < +ka^ i^)^ 



(10) 



where k is related to the confidence level on the pre- 
diction of the limit-state. Any point that belongs to 
this margin corresponds to an uncertain sign of the 
5-function. Thus the scatter of this margin should be 
reduced in order to make the reliability estimation ac- 
curate. This is achieved by adding new points of the 
experimental design in this margin. 

The probability that a point x eVx falls in the mar- 
gin M (with respect to the epistemic uncertainty in the 
prediction) can be expressed in closed-form as fol- 
lows: 



fccTg (a;) — G[x) 



CJQ (x) J 

—kaQ {x) — G {x) 
ag [x) 



(11) 



Two strategies may be used in order to add new 
points: 

• maximise the probability in Eq. ( [TT] ); 

• define a sampling density that concentrates in the 
margin and sample according to this density, as 



proposed in Dubourg et al.| ( |2()T0| ). 

To do so, a sampling density is defined as the prod- 
uct of the probability to fall in the margin of uncer- 



tainty (Eq.( 11)) and some weighting PDF w: 

C{x) (xV{x (^M) w{x) 



The latter equation may be considered as a PDF up 
to an unknown but finite normalizing constant. The 
weighting density w is selected here as a uniform PDF 
on a sufficiently large confidence region of the orig- 
inal PDF. Given that pseudo-PDF, one can generate 
N (say N = 10^) samples by means of any well suited 
Markov chain Monte Carlo (MCMC) simulation tech- 
nique ( [Robert and Casellal |2004| ), such as the slice 
sampling technique ( |Neam2003[ ). 

The N generated candidates are expected to be 
highly concentratred at the modes of C - and thus in 
M. This population of candidates is then reduced to a 
smaller one that uniformally spans the margin M by 
means of the X-means clustering technique ( |EUcan[ 
2003] ). 

Figure [1] illustrates the proposed refinement strat- 
egy applied to a nonlinear limit-state surface. The up- 
per subfigures show the contours of the refinement 
pseudo-PDF C at each refinement step together with 
the candidate population generated by slice sampling 
and its K cluster centers - K ^10 being given. In the 
lower subfigures one can see the real limit-state sur- 
face represented as the dashed black curve, its kriging 
prediction represented as the black line and its associ- 
ated margin of uncertainty M which is bounded below 
by the red line and above by the blue line. Another 
interpretation of these figures is that any point within 
the blue bounded shape is positive with a 95% confi- 
dence and any point inside the red bounded shape is 
negative with the same confidence level. 

In order to stop the iterative refinement algorithm, 
a measure of the scatter of the margin is necessary. 
Introducing the three approximate failure subsets de- 
fined as = {cc : G + i fc erg (cc) < 0}, i = -1, 0, +1, the 
following statement holds: 



The spread of the interval [P(X e F+i), P(X e p-^)] 
is then used in this work as a stopping criteria for 
the proposed refinement procedure. Note that due to 
the inconvenient order of magnitude of these proba- 
bilities, a logiQ-scaling is used in order to assess the 
spread of this interval. Note also that it is neither 
proved or even stated that the real failure probability 
belongs to this interval - due to the empirical assump- 
tions relative to the meta-modeling procedure. How- 
ever the overall strategy has proved to be accurate for 



some applications ( ,Dubourg et al.[ |2010[ [2011) 



(12) 



4 THE PROPOSED ADAPTIVE 

SURROGATE-BASED REDO ALGORITHM 

The proposed strategy to solve the RBDO problem in 
Eq. ([T]) consists in nesting the previously introduced 
kriging surrogate together with the proposed refine- 
ment strategy within a classical (but efficient) nested 
RBDO algorithm. 
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(a) Initialization: 10-point DOE 









(b) Step #3: 40-point DOE 



(c) Step#9: 100-pointDOE 



Figure 1: An illustration of the proposed limit-state surface refinement technique on an example from Waarts ( 2000| l. 



4.1 The augmented reliability space 

The kriging surrogates are built in the so-called aug- 



mented reliability space such as defined in Taflanidis 



and Beck (2009). In this space, the augmented ran- 
dom vector X [6) has a PDF h which accounts for the 
uncertainty in the design choices and the aleatory un- 
certainty in the random vector X. Under such consid- 
erations, h reads as follows: 



h{x) = [ fxix\e)7Ti9) 

JVg 



(14) 



where fx is the PDF of X (6») and tt is the PDF of 9 
that can be assumed uniform on Vg. The augmented 
reliability space is illustrated in Figure [2] in the case 
of a Gaussian variable with a uniformally distributed 
mean. 

The DOE should uniformally span a sufficiently 
large confidence region 'Dx{e) of this augmented PDF 
in order to make the surrogate limit-state surfaces 
accurate anywhere it may be evaluated. More pre- 
cisely, it is expected to be accurate for extreme de- 
sign choices 6 {i.e. located onto the boundaries of 
the optimization space Vq) and extreme values of the 
marginal random vector X (to be able to compute re- 
liability indices as large as e.g. = 8). A confidence 
region is essentially the multivariate extension of the 
univariate concept of confidence interval. Under the 
previous general assumptions, giving a mathematical 




Figure 2: Illustration of the augmented reliability space. 

form to the contour of this region is not as straight- 
forward as for the standard normal space (/3-radius 
hypersphere). However, one may easily build an hy- 
perrectangular approximation by looking at the mar- 



gins of this augmented PDF h as described in Dubourg 



etal. (2010) 



4.2 The adaptive surrogate-based nested RBDO 
algorithm 

The kriging surrogate together with its adaptive re- 
finement procedure is finally plugged into a double- 
loop RBDO algorithm. The outer optimization loop 
is performed by means of the Polak-He optimiza- 
tion algorithm ( Polak[p"997| ). Given an initial design, 
the algorithm proceeds iteratively on i in two steps: 



(i) the direction is determined solving a quasi- 
SQP sub-optimization problem and (ii) the step size 
5^*'' is approximated by the Goldstein- Armijo approx- 
imate line- search rule. The base b in the approximate 
line-search is set to 0.6 and the algorithm usually 
steps forward with < fc < 10 - provided the design 
variables are normalized with respect to their initial 
value. Note that in order to minimize the overall num- 
ber of nested surrogate-based reliability analyses, the 
step size k is saved from one RBDO iteration to the 
other. The reliability and reliability sensitivity anal- 
yses are performed with the subset simulation tech- 
nique ( Au andBeck[ 2001) onto the kriging surrogate. 
The subset simulation technique for reliabihty sensi- 
tivity analysis is detailed in|So ng et al.| ( |2009] ). Sub- 
set simulation is used instead of crude Monte-Carlo 
simulation for the sake of efficiency. Indeed, even if 
the kriging surrogates are faster to evaluate than the 
original limit- state functions, they still involve a non- 
negligible computational effort. Namely, for a 10- 
variable kriging surrogate built from a DOE contain- 
ing a few hundreds points, 10^ runs of the surrogate re- 
quires an average time of 20 seconds for the mean pre- 
diction G and 20 more seconds for the prediction vari- 
ance (7?.. These times were measured on an Intel(R) 
Core(TM)2 Duo CPU (T9600 @ 2.80GHz) running 
with 4GB memory under Ubuntu 10.04 LTS 64bits 
using the built-in vectorization capacity of Matlab 
R2008a. 40 seconds for 10^ complete probabilistic 
predictions is a lot faster than 10^ Finite Element 
runs, though it is still quite time-consuming when in- 
volved in the proposed algorithm. Indeed, the DOE 
enrichment procedure together with the surrogate- 
based reliability analyses (on the three subsets to 
check the surrogates accuracy) are nested within the 
constrained optimization loop. The overall method- 
ology was implemented within the FERUM toolbox 
( [Bourinet et al.[ |2009i ). It makes use of the Mat- 
lab's toolboxes functions quadprog (for the SQP sub- 
optimization problem) and si ice sample (to generate 
samples from the refinement pseudo-PDF C). We give 
a simplified flowchart of the algorithm in Figure [5j 

5 APPLICATION TO THE OPTIMAL DESIGN 
OF AN IMPERFECT SUBMARINE 
PRESSURE HULL 

In this section, the proposed strategy is applied to the 
optimal design of an imperfect submarine pressure 
hull whose mechanical and stochastic models are in- 
spired from pubourg et al.| ( |2008j ) . 



5 . 1 Mechanical model 

The hull is modeled as a single bay of a ring- stiffened 
shell whose axisymmetric section is illustrated in Fig- 
ure |4j This hull is supposed of infinite length and it is 
subject to an external hydrostatic pressure p. Indus- 
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Figure 4: Single bay reference structure - Types and amplitudes 
of imperfections. 

trial structures feature some unavoidable initial shape 
imperfections that are due to the manufacturing pro- 
cess (welding, cold bending, etc. ), to heavy loads 
connected to the hull, etc. Even if these imperfections 
are of moderate amplitudes (thanks to the stringent al- 
lowance for fabrication), they may trigger buckling or 
premature plastic collapse at pressure far below than 
those corresponding to the ideal elastic buckling as- 
sumption. Predicting these collapse pressures is not 
straightforward and is either based on semi-analytical 
formulas (preliminary design) or on high fidelity non- 
linear finite element models (final design). 

The present application is based on the semi- 



analytical formulas presented in Dubourg et al. 
( 2008| ). The code of practice usually demands that 
the stake-holders validate their design with respect 
to three typical buckling patterns known as the over- 
all buckling mode, the interframe buckling mode and 
frame tripping. The two first buckling modes are as- 
sumed to be triggered by two typical shapes of im- 
perfection, respectively: out-of-roundness and out-of- 
straightness. Their amplitudes are denoted by An and 
A„i respectively. For the specific design of interest 
(mean values of Table [T]) n and m values are re- 
spectively set to 2 and 14. These imperfections are 
illustrated in Figure |4] and the associated collapse 
pressures respectively p^pi and p^pi are explicited in 
Dubourg et al.| (|20081). For the last local instability 
though, there is no analytical formula for the collapse 
pressure (at least to the authors' knowledge). How- 



ever, the BS5500 (section 3.6.2.2.2) provides instead 



two arbitrary bounds (using built-in partial safety fac- 



tors) on the stiff eners' proportions ( [Dubourg et al. 
200S1). 



5.2 Stochastic model and failure scenarii 

The stochastic model is given in Table [T| All stochas- 
tic variables are assumed independent. The standard 
deviation of the design variables are kept constant 
along the optimization. For both amplitudes distribu- 
tions, a lognormal distribution with a 50% coefficient 
of variation is assumed. The mean values are deter- 
mined so that the 99.5% quantiles coincides with the 
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Figure 3: The proposed adaptive nested suiTogate-based RBDO approach. 



Variable 


Distribution 


Mean 


St. dev. 


Units 


E 


Lognormal 


200000 


10000 


MPa 


V 


Deterministic 


0.30 






Oy 


Lognormal 


390 


19.5 


MPa 


Ls 


Deterministic 


600 




mm 


R 


Deterministic 


2488 




mm 


e 


Lognormal 


24 


0.72 


mm 




Lognormal 


156 


4.68 


mm 




Lognormal 


10 


0.30 


mm 


Wf 


Lognormal 


120 


3.60 


mm 


e/ 


Lognormal 


24 


0.72 


mm 




Lognormal 


5i?/3000 


0.5 X 5i?/3000 


mm 




Lognormal 


L^/300 


0.5 X L,/300 


mm 


Pq 


Deterministic 


2 




MPa 



Table 1: Stochastic model for the single bay reference structure. 



maximum allowable values An2 = 5i?/1000 (BS5500I 
and An 14 = is/100 (choice of the authors). Due to 
the lognormal assumption this is achieved for a mean 
value equal to about the third of the targeted 99.5% 
quantile. 

Two failure scenarii are considered in this paper. 
They are defined with respect to an arbitrary reference 
diving pressure po ^2 MPa. The first failure scenario 
considers overall buckling and is thus described by 
the following limit- state function: 



5i(X(0))-log 



Pnpl 
PO 



(15) 



and the second failure scenario considers interframe 
buckling so that its limit- state function reads: 



52(X(0))=log 



Pm pi 
PO 



(16) 



The logarithmic scaling of the two previous limit- 
states is chosen so as to simplify the meta-modeling 
task. Finally, the structure is considered to fail when 
either one or both of the two previous scenarii oc- 
curs (serial system combination), so that the limit- 
state function reads: 



g{X{e)) = mm{gi{X{e)), g2{X{e))} 



(17) 



Note that in the proposed strategy a surrogate is built 
for each limit-state, and the final surrogate limit-state 
is evaluated as the minimum of these two. 

5.3 RBDO problem formulation 

The reliability-based optimization is carried out with 
respect to the mean dimensions of the stiffener: ^e, 
Me,„, Me/, i^wf The interframe distance Ls and 
the mean radius R are supposed to be imposed by 
other constraints (nuclear reactor size, etc. ). The opti- 
mal mean design should minimize the expected ''hull 
weight" I ''water displacement" ratio which is approx- 
imated as follows (using the mean values of all the 
random variates): 



E[c(X(0))] w c(E [X(6»)]) 

PsteelVsteel(E[X(g)]) 
Psea water Vsea water [-X^(0)]) 



(18) 



where Vsteei and Vsea water Stand for the hull volume and 
the volume of displaced water, and psteei = 7850 kg/m^ 
and psea water = 1026 kg/m'^ are the densities of steel and 
sea water respectively. Finally, in addition to the reli- 
ability constraint I3{d) > Po, we consider the two de- 



Variable 


Initial 


DDO 


RBDO (Present application) 


RBDO(Dubourgetal. 


2008 , Table 6) 


Units 


= 3.00 


= 4.50 


= 6.00 


(Jq = 3.00 


f^o = 4.50 


(3q = 6.00 




24.00 


16.90 


21.46 


26.10 


32.01 


*20.67 


* 24.47 


*29.82 


mm 




156.00 


160.27 


171.25 


178.35 


189.22 


*161.99 


*197.71 


*231.82 


mm 


Me™ 


10.00 


7.16 


10.15 


10.46 


19.66 


*9.40 


*7.96 


*9.67 


mm 


f 


120.00 


81.89 


89.46 


95.14 


103.36 


* 110.58 


*99.08 


* 116.98 


mm 


^J■ef 


24.00 


16.76 


21.28 


26.04 


31.85 


*22.18 


*24.14 


*27.01 


mm 


C 


18.86 


12.75 


16.57 


20.01 


26.07 


16.41 


18.64 


23.14 


% 




4.04 


-0.44 


3.04 


4.48 


6.01 


2.98 


4.47 


5.96 




C.o.V. 


3.34 


0.02 


2.52 


3.82 


5.40 


2.50 


3.86 


5.30 


% 


51 -calls 




183 


200 


300 


350 


**3257 


**9494 


**6580 




g2-calls 




183 


200 


300 


350 





* The design variables were considered as deterministic. ** The authors provide only the total number of calls. 
*** Reliability indices checked on the real limit-state functions using subset simulation with N = 10^ runs/step. 
Table 2: Optimal designs for the single bay reference structure. 




Iterations Iterations Iterations Iterations 



Figure 5: Convergence of the proposed algorithm for the minimum reliability constraint /3o = 6.00. 



terministic safety constraints from the BS5500 given 



in 



Dubourgetal. (2008) 



5.4 Results 

First, the starting design for the RBDO problem is de- 
termined by solving the equivalent deterministic de- 
sign optimization (DDO) problem. To do so, all the 
random variates are replaced by their mean value and 
the probabilistic limit-states are used as deterministic 
constraints. The deterministic optimal design is e = 
16.90 mm, hu, = 160.27 mm, = 179.86 mm, Wf = 81.89 
and = 16.76. This design lies on the boundaries of 
the two limit- states gi = and 52 = and is conse- 
quently rather sensitive to uncertainty as it is associ- 
ated with a generalized reliability index of /3 w -0.47 
{i.e. a failure probability greater than 50%). 

Then, the proposed strategy is applied to solve the 
reliability -based design optimization starting from the 
deterministic optimal design and for three different 
minimum reliability constraints: = {3.00, 4.50, 6.00}. 
The algorithm uses a 50-point initial DOE and it is se- 
quentially enriched with 50 new points if required - 
i.e. if the logiQ-scaling of the spread between V{X e 
F+i) and V{X e p-^) is greater than ep^g = 5 x lO'^. 
The results are given in Table [2] together with the ini- 
tial design (which corresponds to the current prac- 
tice), the deterministic optimal design and the FORM- 



based optimal designs from Dubourg et al. (2008 1. 

For each design, we provide the associated ob- 
jective function value, the generalized reliability in- 



dex estimated by subset simulation (using the real 
limit-state function) with its coefficient of variation 
(C.O.V.), and the overall number of calls to the limit- 
state functions. None of the local instability safety 
constraints prescribed in the |BS5500| are saturated at 
these optimized designs. Note that even if the number 
of calls to the limit-state functions increases with the 
minimum reliability target, convergence is achieved 
using a few hundreds evaluations of each only com- 
pared to a few thousands for the FORM-based double- 
loop approach proposed in 'Dubourg et al.|([2008). The 



reliability-based optimal design provided in this pa- 
per are significantly less optimal than the one given 
in Dubourg et al. ( 2008j ). This may be explained as 
follows: 

• the design variables were assumed deterministic 
in Dubourg et al. ( |2008) whereas they are con- 
sidered as random with a 3% C.o.V. (w.r.t. the 
initial mean design of Table [T]) in the present ap- 
plication; 

• the FORM assumptions do not seem to be con- 
servative in this case since the reliability indices 
estimated by subset simulation are lower than the 
minimum reliability constraints. 

The convergence of the proposed algorithm is de- 
picted in Figure |5] for the minimum reliability con- 
straint /3o = 6.00. It can be seen from the extreme- 
right subfigure that the number of limit- state function 
calls increases as the optimization algorithm moves to 



higher reliabihty levels. It makes sense as the lower 
the failure probability is, the narrower the margin 
of uncertainty has to be. The extreme-left subfigure 
shows the convergence of the design variables nor- 
malized with respect to their starting deterministic op- 
timal values. It means that the product of the optimal 
characteristic mean design (from DDO) by these op- 
timal normalized values would ensure a ^o-reliabiUty 
optimal design. They could thus be interpreted as op- 
timal partial safety factors, though they were cal- 
ibrated with respect to a single reference structure 
{R = 2880 mm, Lg = 600 mm) as code calibration usu- 
ally requires to consider several reference structures. 

6 CONCLUSIONS 

The proposed surrogate-based approach allows to 
save a significant number of simulation runs without 
using the first-order reliability theory assumptions. 
It is important to note that its efficiency mainly re- 
lies upon the fact that the kriging surrogates are built 
in the so-called augmented reliability space and that 
they can thus be used from one nested reliability anal- 
ysis to the other. The results provided in this paper 
provides some reference results for a future applica- 
tion to the design of imperfect submarine pressure 
hulls that will involve a higher fidelity mechanical 
model, namely a nonlinear shell finite element model. 
Indeed, it is of the authors' belief that reliability-based 
designs have to be based on the best up-to-date nu- 
merical methods, for meaningful and trustworthy re- 
sults, and without depending too much on the degree 
of conservatism inherent to closed-form or approxi- 
mate solutions such as the one used in the present pa- 
per. 
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